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We consider the problem of two-level system dynamics induced by the time-dependent field B — 
{a{t) cos (jjt,a{t) sin (jjt,(jJo}, with a{t) oc cn{i't,k). The problem is exactly analytically solvable and 
we propose the scheme for constructing the solutions. For all field configurations the resonance 
conditions are discussed. The explicit solutions for TV = 1, 2 we obtained coincide at tj = in 
the proper parameter domain with predictions of the rotating wave approximation and agree nicely 
with numerical calculations beyond it. 
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I. INTRODUCTION 

Dynamics of model two-level systems in the external 
fields is of fundamental interest, as this problem appears 
in various physical systems. They are, for instance, spin 
1/2 in the magnetic field two-level atom interacting 
with classical field [H, driven artificial two-level systems 
in superconducting Josephson devices Q. As well this 
problem is significant for physical realization of quantum 
gates [1] . But still there are few time-dependent field con- 
figurations for which the problem is exactly analytically 
solved 0-01) so in many cases the solutions are given us- 
ing perturbation theory or computed numerically (8l-[To| . 

One of the physically significant fields is linearly po- 
larized (LP) harmonic wave. In particular, it is used in 
magnetic and optical resonance experiments [H, 0] . There 
is a number of approximate approaches to this problem. 
In particular, when one assumes the amplitude of the os- 
cillating field to be small in comparison with the static 
field along the third axis and the frequencies of these 
fields are near-resonant. In this case the rotating- wave 
approximation (RWA) is valid and one can easily get the 
solutions known as Rabi oscillations (l| . The understand- 
ing of magnetic resonance, meaning the resonance con- 
ditions, is based on these approximative treatment, and, 
for example, the resonance frequency is given by series in 
the ainplitude (the first correction is called Bloch-Siegert 
shift) [nl, [13] ■ On the other hand, not so long ago this 
problem was shown to be exactly analytically solvable 
in terms of Heun functions ?| . But as the apparatus of 
the latter is not developed enough, it is quite difficult to 
extract physical information from the solutions. 

In this Letter we consider the generalization of the 
Rabi field which at different values of parameters corre- 
sponds to several physically significant situations: non- 
linearly modulated Rabi field, N-soliton pulse and the lin- 
early polarized field, which approximates harmonic one. 
We succeed to obtain explicit analytical expressions de- 



scribing dynamics of a two-level system in the proposed 
time-dependent field and, consequently, to get conditions 
of resonance in it. The emphasis should be laid on the 
third field configuration, as it refiects the idea to replace 
the cosine potential with another one resembling it with 
arbitrary accuracy but corresponding to the exactly solv- 
able problem in terms of well-developed functions. The 
exact results we obtained coincide at small amplitude of 
the field and near-resonant frequencies with predictions 
under the RWA, while beyond it they are in agreement 
with numerical calculations of the problem with the lin- 
early polarized harmonic field. The resonance conditions 
in the analytic form for this field, in particular, reproduce 
at small field amplitude the Bloch-Siegert shift. 



II. FORMULATION OF THE PROBLEM 

We consider the problem of a two-level system dynam- 
ics in the external field 



B(i) ~ {2a{t) cosujt,2a{t) smujt,uJo}, 



(1) 



where a{t) — Nkvcx\{vt,k), N ^ "L. Here Civ{vt,k) de- 
notes the Jacobi (cnoidal) elliptic function of modu- 
lus /c, < A: < 1 [l^. Let ut analyze briefly the 
range of physical situations given by the field B(i) 
For the arbitrary parameters it can be understood as 
nonlinearly modulated Rabi field, meaning the nonlin- 
ear function a{t) rather than constant amplitude a in 
B_R(t) = {2acoswt, 2asinwt, Wo} [l|. In the case k = 1 
the amplitude is a{t, k = 1) = Ni/sechi/t, so the field is 
rotating with the iV-soliton envelope 3] . The third situ- 
ation occurs under uj — : BLp(t) = {2a{t), 0, luq}- If the 
limit iV — T' oo, A: — ^ {Nk = const) is taken, this field 
is exactly the linearly polarized harmonic wave. There is 
also an approximate limit, valid for appropriately small 
k, when cn(j/t, fc) « cos( 2K(k) ^(^^^ K{k) is the com- 
plete elliptic integral of the first kind. As the period of 
the cnoid is AK{k), due to definition of the frequency 

we obtain = vk /\2.K{k)] as an "effective frequency" 

of the cnoidal field. Therefore, finding the dynamics in 
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small enough k gives solution to the problem with the 
field {cosfii, 0, wq}- 

The Hamiltonian of our problem has the form: 



(2) 



with a — {a-i} - Pauli matrices. Let |±) be the eigenstates 
of operator with eigenvalues ±1, respectively. The 
evolution of a generic state of a two-level system |^'(t)) = 
where \C+{t)\^ + \C-{t)\^ = 1, is 
determined by the Schrodinger equation 



(3) 



Considering the form of Hamiltonian ([2]) for the given 
choice of fields, we have the system of equations for am- 
plitudes C±{t): 
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We will choose the initial conditions as 

C+(0) = 1, C_(0)=0, 

meaning that the dynamics starts from the 
state. 



(4) 

(5) 
-eigen- 



III. SCHEME FOR CONSTRUCTING THE 
SOLUTIONS 



Introducing new variables C±{t) = e^^'^°*/'^C±{t) the 



system of Eqs. ^ can be reduced to independent equa- 
tions of the second order. In terms of the dimensionless 
variables they take the following form 



snr dnr 



±iA 



dC± 



-iV^/c^pnV C± = 0, (6) 



where snr, cnr, dnr are the Jacobi elliptic functions of 
modulus k, T = vt and A = (ljq — uj)/v is dimensionless 
detuning. After the change of variable s = sn(T, k) the 
equations © are 



d^C± 
ds^ 



1 - fc2s2 



± 



dC± N^k^ 



ds 1 — fc^s^ 



C± = 0, (7) 



thus they can be classified as the equations of the Fuchs 
type J^j on the Riemann surface corresponding to the 
algebraic curve w"^ = {1 — s^)(l — fc^s^) 13]. 

Two obvious limiting cases (A = and fc = 1) lead to 
simplification of the Eqs. (O to the hypergeometric type 
which could be treated easily and are physically mean- 
ingful. In the case A = 0, which means the resonance 
between frequencies of constant and oscillating fields, the 
Eq. ([7]) is simplified to equation of hypergeometric type, 
with solutions: 



Tjv(-fcs), C- = - k^s^\JN-i{~ks) (8) 



where Tn{x) and Un{x) are Chebyshev polynomials of 
the first and second kinds, respectively. The amplitudes 
satisfy proper initial condition automatically. Other sit- 
uation can be obtained setting fc = 1, so that initial field 
is reduced to the rotating field with N-soliton envelope 
a{t,k = 1). The equations are also hypergeometric and 
the solution is 



(9) 



where P^'^^ is the Jacobi polynomial with complex pa- 
rameters a = —(1 — «A)/2, b — —(1 -I- iA)/2. Here it is 
important to mention that the initial condition should be 
substituted by C+(— oo) = 1, C-(— oo) = 0, as the field 
is non-periodic and its action starts aX t ~ —oo. 

Let us seek for solution of the equation ([7]) when " plus" 
is chosen in the form 



/+ = exp 



7^o(s')+^ 



7^l(s' 



ds' 



(10) 



where 7^o(s) and TZi{s) are rational functions in s. As a 
result, the initial equation ([7]) comes up to the system of 
Riccati-type equations 



n' 



7^i + 7^l(27^o - 



(ATCi-TCj) 



A7^o = 0. 



O(lla) 
(lib) 



Obviously, TZo{s) in terms of the amplitude modulus 
has the form 



Since solutions in both limiting cases imply that |C-| 
a polynomial function, we will try in the form 



(12) 



N 



where V2N is a polynomial of the 2A^-th order with roots 
Ci. Then due to TIq{s) can be easily calculated. The 
solution of (|llb[) with known TZo{s) is given by 



A I e 
'0 



The equation (|llap serves as the consistency condition 
and defines the roots e^. With known TZo{s) and TZi{s) 
one can obtain the expression for the /+ due to (|10l) . 
Analogically the solution /_ of ([7]) standing for the "mi- 
nus" sign can be found. Then one needs to chose two 
linearly independent solutions of ([7]) and create linear 
combinations satisfying initial conditions 

Using this scheme we can find the solutions for = 1 



C+ = £i^s2 + e2 c-"^' +82^-8^ -ej e-*'^^ (13a) 
C_ = -i£2a/s^ + e? e'^P' +ieiJ-s^-el e"f\ (13b) 
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where e\ 2 = 



A2 - 1 ± 



/(2fc2 



£j = V(-l)*+^ef/(e? - e^), i = 1,2. Here the time- 
dependent phase 

= A [r - (1 + e,2)(eD-^ H (-(e^)-!; am(r, fc), k)] 

is defined in terms of the incomplete elliptic integral of 
the third kind 

n(n;a m(T, fc),fc) = /^[(l 

nx^)\J (1 — a;2)(l — fc^a;^)]"^^^: and am(T, fc) denotes the 
Jacobi amplitude function [Tsj . 

Increase of (physically it corresponds to the increase 
of the field amplitude) results in the complication of ex- 
pressions for the amplitudes but doesn't change its struc- 
ture. Explicit solutions for iV = 2 are the following 



C+(t) - eav/(-'S' + ei)(s2-e2)e'^(^-'=^) 

+ e6v/(-s' + e3)(-s2 + e^y^^'^-''^\ (14a) 



C_(t) - iebV(-s' + ei)(s2-e2)e-*'^('=i^^^' 

- i£aV(--'^' + e3)(-s2 + e4)e-*'^(^^'^^). (14b) 



Here £a = ^-6162/(6364 - 6162), = ^6364/(6364 

6162), 

e, = [-4(A2 - 9) (-1)^+1n/2A±] /{36P), ^ = 1,4, 



9 (9 ± + 54 (1 - 2fc2) A2 - 7A4 



1/2 



where "plus" stands for 61^2 and "minus" for 63^4, 
d = [(A2 - 9) (A2 - 1) + 12(Afc)2] ' + (A2 - 9 + 48fc2)' 
x(4A/c)2 and the phase 

ip{ei,en) = A [2r/3 - 7;H; -t- 7„H„] , 



7i(n) = (ei(„) - l)(-4 + e;(„) + 3e„(i))/ [36;(„)(6i - 62)] , 
where {l,n) stands for the pairs (1,2) and (3,4), and 
=n(e,-^am(T),fc) . 



IV. RESONANCE CONDITIONS 

If the amplitudes are known, one can easily calculate 
the expectation values of spin momentum projections 
which are defined as 8^(1) = {9{t)\ai\^{t)) . They stand 
for the Bloch vector vector S{t), whose evolution takes 
place on sphere. Its third component S3 {t) , which is use- 
ful for physical interpretations, can be found in terms of 
the amplitudes by the expression 



S,{t) = \C+{t)\^~\C^{t)\^^ 



(15) 



where the initial conditions ([5]) are taken into account. 

The problem of two-level system dynamics in the Rabi 
field Bfl, which approximates the linearly polarized har- 
monic field when the amplitude of latter is small and the 
frequencies are near- resonant (RWA) , can be easily solved 



[T|. The third component of the Bloch vector oscillates 
(Rabi oscillations) due to 



4a^ cosnBt]/^% 



(16) 



where Qn — vA^-|-4a2 is the Rabi frequency and 
A = — is the detuning between the rotating and 
the constant field frequencies. 

The term resonance we will refer to the situation when 
under some special relations involving characteristic fre- 
quencies, system evolves from |-|-) to |— ) state, thus min- 
imum of 5'3(t = tq) is — 1. In the Rabi field B/} the only 
resonance condition is w = wq, then the system performs 
oscillations with the period Tr = n/a between states. 
The situation when the characteristic period of the dy- 
namics is Tji we will designate as the Rabi resonance. 

In the nonlinearly modulated Rabi field uj = ujq is not 
enough to satisfy S'3(to) = —1. Using the expficit expres- 
sions for the amplitudes given by ([8|) we obtain 

Ssir) = 2T7v(-fcsnT)2 - 1. 

Thus TAr(— /cres)^ = is the additional resonance con- 
dition (Fig. [T]). For the situation of the N-soliton pulse 
(k — 1) the only resonance condition is v = loq. 



(b) 




FIG. 1: Resonance in the nonlinearly modulated Rabi field, 
(a) iV = 2, kres = 0.7, (b) iV = 3, kr^s = 0.89. 

In the case of two-level system driven by B^p (r) , S3 (r) 
can be calculated for iV = 1, 2 using explicit expressions 
(1121) -(dH). Our analytic results are in good agreement 
with numerical calculations for the linearly polarized har- 
monic wave in wide range of parameters (see Fig. [5]) . As 




FIG. 2: Comparison of the exact result for N = 2 (solid 
line) in the field BLp(r) with the numerical integration (cir- 
cles) of the problem with the linearly polarized harmonic field 
{2acosr,0,tJo}, a = 2k. Here k = 0.25, a = 0.5, (a) A = 0.4, 
(b) A = 12. 
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the field Blp(t) has two characteristic frequencies: ujq 
and ri, the expected resonance condition is 



(17) 



2K{k) 

For small k the last expression can be expanded as series 

5fc4 



64 



(18) 



For = 1, fc plays the role of the constant amplitude 
a in the corresponding B^j, thus lllSp coincides with ex- 
pression for the Bloch-Siegert shift [lllfl^. 



(a) 



(b) 





(c) 



(A) 





Rabi resonance occurs under 



fc 2^1(1/2, 1/2; l;fc2) 



mN 



(19) 



The solutions presented here under resonance condi- 
tions = llIq and (jl9p are consistent with the results of 
the numerical integration of the problem with linearly 
polarized harmonic wave of the same amplitude and fre- 
quency. As well, they resemble dynamics predicted under 
the RWA (Fig.[3|). If one is interested to study dynamics 
of the two-level system induced by the linearly polar- 
ized harmonic field in the extended range of amplitudes 
of time-dependent field, it is necessary to calculate solu- 
tions with higher N due to the proposed scheme. As the 
amplitude is TVfc, the increase of N also allows to decrease 
k making the approximation better. 

If under a two- level system the spin 1/2 system is re- 
garded, the results obtained here can be easily general- 
ized to the spin j system (equivalently, to 2j 1 non- 
interacting spins 1/2). The amplitudes C±{t) with re- 
spect to the initial conditions ([5]) form the matrix of evo- 
lution u{t) S SU(2), with the elements a — C+, (3 = 
~C1. With known u{t) one can get the solution for the 
spin j system due to the well-known procedure (see for 
detail m)- 



V. CONCLUSIONS 



FIG. 3: Comparison of the exact result for N = 2 (solid line) 
in the field B_Lp(r) with the numerical integration (circles) of 
the problem with the linearly polarized field {2a cos r, 0, ljo}, 
a — 2k. Here = ujq and k are chosen to be solutions of 
the k 2-Fi(l/2, l/2;l;fc^) = (2m)"\ therefore, the Rabi os- 
cillations (dashed line) with period Tr = 2n/a under small 
k are expected to occur, (a) k = 0.1245, (b) k — 0.1655, (c) 
k = 0.2461, (d) k = 0.4701. 

Our results are expected to be consistent with Rabi 
theory when RWA works. In the terms used here it cor- 
responds to the approximate resonance condition v = ujq 
(A = 1). For A*" = 1 under this relation the expression 
for the S3{t) can be simplified 



1 , / dnr — zfcsnrcnT 
fcr + - In 

2 ydnr + zfcsnTcnr 



In the corresponding Rabi field Bp dynamics is given 
by S3{t) — cos/cT, as a = fc. Therefore, simultaneous 
satisfaction of sn(To,fc) = and cos(fcTo) = 1 (tq = 
2mK{k), m E Z and tq — 7r(2n + l)/fc, n £ Z, re- 
spectively) gives values of k corresponding to Rabi reso- 
nance. The ratio m/(2n + 1) have to be an integer, we 
set n — for brevity. Using the known relation between 
the incomplete elliptic integral of the first kind K{k) and 
the hypergeometric function jl3l | , the resonance condition 
for TV = 1 takes the form k 2^1(1/2, 1/2; 1; fc^) = l/m. 
It can be generalized for the arbitrary N, meaning that 



In this work we have considered the problem of two- 
level system dynamics induced by the external time- 
dependent field, expressed in terms of the Jacobi elliptic 
function (cnoid) . We presented the scheme for construct- 
ing the exact solutions of the Schrodinger equation with 
such a field and used it to calculate explicit expressions 
for the amplitudes. At different parameters the chosen 
field corresponds to three physical situations: a nonlin- 
early modulated Rabi (rotating) field, with fc = 1 it is 
the N-soliton pulse and with lu — and small enough 
k it resembles linearly polarized harmonic wave. For all 
these field configurations we discussed the resonance con- 
ditions. In particular, in the case of linearly polarized 
field the resonance condition in the analytic form repro- 
duce at small field amplitude the Bloch-Siegert shift. The 
explicit solutions for = 1, 2 with a; = are shown 
to coincide at small field amplitudes and near resonance 
with predictions of the rotating wave approximation and 
to be in good agreement with numerical calculations in 
the wide range of field parameters beyond it. 
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